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Q\ ; I. INTRODUCTION 

ON ' 

Recently [Q-p| the behavior of the spatial string ten- 
sion was studied in (d+ l)-dimensional (d — 2,3) SU(2) 
and SU(3) gauge theories. The spatial string tension a 
at high temperature in (3 + l)-dimensional SU(N) gauge 
theory was rigorously proved H to be non-vanishing at 
finite lattice spacing. The spatial string tension, as it was 
pointed out by , was not related to the confining prop- 
erties of a physical potential in the (3 + l)-dimensional 
theory. The reason is that under Z (N) - transforma- 
| tion topologically trivial Wilsons loops remain invariant; 

on the contrary, topologically non-trivial loops such as 
| Polyakov lines transform as (x) — > zCl (x) . Therefore 
the behavior of topologically trivial Wilsons loops cannot 
be considered as a confinement criterion [Q. In partic- 
ular, the expectations of a large space-like Wilsons loop 
may show area law behavior without static quarks being 
confined. 

Despite the fact that the study of space-like Wilsons 
loops behavior at finite temperatures does not give 
■ straightforward information about critical phenomena in 
LGT, it helps to understand better non-perturbative ef- 
• i-h . fects that manifest themselves in correlation functions for 
the spatial components of gauge fields. 

The remarkable feature of the spatial string tension is 
that it is scaling of j| and thus ^ is non- vanishing in 
the continuum limit. The calculations of an average value 
of a time- and space-like Wilsons loop at g 2 ~ gcritical i- e - 
at fixed cut-off have been performed for a broad temper- 
ature interval (T was varied by varying N T ) for SU (2) ||] 
and for SU (3) [Q gauge groups. It was shown that the 
spatial string tension remained temperature independent 
up to T c and than was rising rapidly, unlike the temporal 
one that decreased with temperature above T c . Similar 
behavior has been found in lower dimensions and also in 
Z{2) gauge theory 

The main features of high temperature behavior of 
such observables as the heavy quark potential and spatial 
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string tension can be understood in terms of the struc- 
ture of the effective, three-dimensional theory which was 
obtained from dimensional reduction at high tempera- 
ture by means of perturbation theory . The basic sug- 
gestion is that at high temperatures temporal dimension 
becomes arbitrary small and degrees of freedom in that 

direction are frozen, therefore X^tllo Ex S ~ ITT Ex S 
so for three dimensional couplings one can write 

with 

g 2 m =g 2 (T)T. (1.2) 

As it was established in j|, the Higgs part of such 
an effective theory does not contribute substantially in 
the spatial string tension, leading to the simple relation 
between 03 and a. One of the remarkable results of 0] 
can be given as 

V^~cN-gf 3 y, N = 2;3. (1.3) 

2 

It means that y/^/T c cN (3 jT" = cNg 2 rlt is in agree- 
ment with the MC experiment, which shows that the 
scaling violation of the ratio y/a/T c is small enough Q. 

Though the results obtained in a MC simulation have 
already answered many crucial questions, we hope that 
an attempt of an analytic study presented here will be 
useful for a more detailed understanding of a scaling phe- 
nomenon. The present paper is organized as follows. In 
Sect. 2 we discuss the main suggestions, that have been 
made in the model and compute the average value of a 
spatial Wilson loop in given approximations. In Sect. 3 
we discuss the result, obtained for spatial string tension. 
In Section 4 we make an attempt to give a more compre- 
hensive discussion of the spherical model approximation 
accuracy and list some standard and nonstandard exam- 
ples of model applications. 



II. MODEL 

To compute an average value of the space Wilson loop 
WnxL((3a) = (Wijxi), we shall use the asymmetrical 
lattice — = £ ^ 1 so for the action one can write 

S = Se + Sh', 

JV T -l 1 

Sh = -0* J2 ^ReS P D m „ (2.1) 

t— x,t,niy^n 
JV T -1 



S E = ~PtJ2Y, ^ReSpn „; (2.2) 

with 



N 

t—0 x,n 
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2N 2N£ 

Pa = Pnm = ~T7'> Pt — Pon — — 5~ • (2-3) 

git 9r 



and 



□mv = ^ (s) ^ + m) (a + E# (a) (2.4) 



Lattice spaces ao and a a can be made arbitrary small at 
any fixed £ and no special assumptions are made about 
g% and g%. except that their values are wholly determined 
by certain renormalisation group equations at any given 
a T , a a , N T , and with N% = N1N2N3. With decreasing 
£ the electric part (Se) of the action becomes negligibly 
small in comparison with the magnetic one {Sh), then in 
£ << f limit (which in a way is opposite to the Hamil- 
tonian limit: £ >> 1) it may be ignored. The magnetic 
part of the action in this case is split into a set of indepen- 
dent time slices. In other words, the temporal degrees of 
freedom are frozen in the limit £ << 1 (which is natural 
at high temperatures 0), and we have 



Sh = Pa -^ReSpD mn 

~Pa-N T ^RcS P D m „. (2.5) 



The Wilson loop is placed in one of the slices t = to and 
is not affected by any other slice, therefore this specific 
part (t = t ) of the action works as an effective action 
while calculating the average value of Wilson loop. To 
put it differently 



Z = J] Sp t exp i A J2 RcS P D »»> )\= Z (*o) 



(2.6) 

and Pa = rj^ can be regarded as the effective coupling. 
Therefore 

W RxL (p a ) = (2.7) 
'' -Sp to I W flxL exp i Y ReSpD„, 



Z(t ) > t0 — 1 ) gU a 



A. Approximation SU(N) ~ Z(N) 

It is commonly believed that the Z(N) degrees of 
freedom arc responsible for many important aspects of 
SU(N) gluodynamics phase structure. The lattice gauge 
theory in the vicinity of phase transition is widely known 
to show large degree of universality 0. The main evi- 
dence in favor of universality is given by the Wilson loop 
behavior, whose functional form does not depend on the 
choice of a gauge group or on the specifics of ultraviolet 
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behavior of the model, showing rather simple dependance 
on space dimension. Universality arguments place the fi- 
nite - temperature SU (N) - gauge theory in the univer- 
sality class of globally Z(N) invariant systems with short 
- range interactions. Hence it is convenient to study such 
universal infrared behavior in SU (N) ^ Z (N) approxi- 
mation: 

U x -f, ~ z z;M = exp | 2m ^^ | ■ q x . lt = Q ) ... ) N-l. 

(2.8) 

The MC experiment as well as model calculations Jh]] 
demonstrate that SU (2) and Z (2) spectra show remark- 
able agreement not only between the pattern of the 
states, but also between the values of the masses (except 
for the lowest state) (see also || and ||). Such mod- 
els with discrete gauge groups are easier to handle. In 
particular, the method of duality transformations is just 
the one elaborated well enough for the systems with dis- 
crete symmetry ( and references there) These models 
with Z(N) gauge symmetry are also known to provide a 
transparent realization of 't Hooft algebra of order and 
disorder operators. Moreover, there is evidence jl3| that 
the effect of the quantum fluctuations near to the Z(N) 
configurations on the symmetrical lattice leads only to a 
finite renormalisation of the coupling constant: 

N 2 - 1 

Pold — * 0new = Pold ^ • (2.9) 

Since the additional term in ( |2.9| ) depends neither on a a . T 
nor on (3 we may hope that on the asymmetrical lattice 
the effect of the quantum fluctuations will also lead again 
to an insignificant change of the coupling constant. We 
consider it to be still another justification of the chosen 
approach. Of course, nobody expects full coincidence of 
the results for SU (N) and Z (N) . 

One of the main advantages of the SU (N) w Z (TV) 
approximation is that one can apply it to Z (N) du- 
ality transformations, which relate Z (TV) gluodynamics 
to Ising (N = 2 ) and Potts (N = 3 ) models in 3 - 
dimensional space . 

B. Duality transformations 

As it is well known [O, in the case of Z(2) gauge group 
the Wilson loop R x L average value placed at x\ = 
can be calculated on the dual lattice 



W RxL = 




z s e Z{2). 



(2.10) 

Summation 53hxl IS done over all dual Ising spins placed 
at x\ = inside R x L. In other words, ferromagnetic 
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links dual to plaquettes of R x L transform into antifer- 
romagnetic links of the same strength JlJ]. In the Z(3) 
gauge theory where the coupling (3' is multiplied by e~s~ l 
at links dual to plaquettes placed inside R x L at x% = 
we have 



W, 



B,xL 



exp < - (l - e^ 4 ) Re ^ /3^%+ r , 



RxL 



(2.11) 



The duality transformations on the asymmetric lattice 
(ai ^ fl2 ^ 03 and Pi 7^ P'2 7^ P3) relate dual links 
l n to plaquettes P m k of the original lattice (n =/= k =/= m) 
p0| . The couplings (3' n on the dual asymmetric lattice are 
related to the corresponding ones {(3km) on the original 
lattice in the form 

Li. ' 7 m 2 ' (2 - 12) 

In the specific case of a n = a, {(3' n — (3' a ; (3km = At ) 

W; ■*<<!: (2 ' 13) 



C. Spherical model approximation 

To compute the partition function^] 



Z = eXP \ Re E ^n^x+ri + E f 

{z} I a,n a: ) 



(2.14) 



on the dual lattice we use the well-known spherical model 
jl7j (see e.g. |Q, JljJ and references therein). The crucial 
point of this model lies in replacing the precise condition 
I z x 1 2 = 1 by a less burdening " averaged" condition 



iV3 



or 



n^(i^i 2 -i)-^(E^ 3 



uj 2 



\z x \ 2 \\, (2.16) 



where the constant c is chosen to the right of all sin- 
gularities of the integrand to ensure the correctness of 



1 Sources r) x are introduced for convenience. 
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interchanging the integration over z x and ds. Now we 
can change — > dz x , so the partition function 



(2.14) may be rewritten in the following form: 

Z - / —<"■"' x (2.17) 



2m 



/ y\ dz x exp - z x Ag z x > + 



where 



3 



" l °x+n 



+ n 



(2.18) 



with 



\ 



(2.19) 



The integration over can be now carried out 
1 / ^-l 



dsenplsN* - -Indet A- r] x [A% ) r) x > } = 



dse^p{N^(s)} 



(2.20) 



To compute the Wilson loop average value, in the par- 
tition function Ai = sS? - J2 n =i P'rJ>3+n sn °uld be 
changed into At + (A RxL ) x s , where 



X 3 + - 



L 

aw 

2 J 



(2.21) 



After the Fourier transformations 
^^exp(iff-i^^).4=A( W ^) (2.22) 



we have 



with 



A(^)^E^ C0S 



(2.23) 



(2.24) 



n=l 



and 



G 



. (R+l)(<p 2 -<p' 2 ) . (Wfe^) 

A(^;^) = 2 . ^ ■ (2.25) 

The integration over s can be done by the steepest 
descent method, the saddle point s being defined by the 
condition [^^(s)] s=So = which can be rewritten as 

1 = 1 r i . (2.26) 

In the symmetric case the equation ( ^.26| ) has a simple 
solution in the critical point vicinity 

so « 3/3' + 27r 2 p'((3' c - I3') 2 (2.27) 

with 

3 

:-! 



m=y / So - 3/3' w ny/WiP'c - PWc - P)- (2-28) 

The number of sites spanned by the Wilson loop is 
I/N2N3 times smaller than the whole volume, so their 
contribution doesn't influence the saddle point position. 

To clarify the e ngin e of spherical model approach we 
would note that ([2. la) fixing the compactness condition 



\z\ = 1 with shown approximation brings (through (2.16) 
) an additional 'mass' term into action: 



3 

m 2 = ao GO-XX. (2-29) 
n=l 



The effecti ve m ass ( [2.28 ) is defined by the saddle point 



condition (2.26) and plays the role of the screening mass 
(expressed in lattice units) in the correlation functions. 

As it is well known, theories with the local gauge sym- 
metry are described in terms of nonlocal order parameter. 
Thus the partition function also may have no singular- 
ities in the thermodynamic limit, the quantities which 
determine the nonlocal order parameters may have sin- 
gularities as a result of increasing their size to infinity. 
This is strongly suspected to occur in lattice gauge the- 
ories for Wilson-loop order parameters, and poses an 
obstacle to the strong-coupling expansion pl| . There- 
fore, it seems quite useful to study differences in behavior 
of a 'tiny' (one-plaquettes) and 'large' (whole N 2 x N 3 
plain) Wilson loops. In particular for the 'small' one 
(1 << R « N 2 ; 1 << L « N 3 ) it would be enough 

to take into account the term ^^'f"^ cos ipi in the first 
order 

RL^frcosvx (2 30) 



~"a73 2^ 



v 
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therefore from 

Wr.xl = cxp | -RLa 2 a 3 a {1) } , 
we have for the string tension 



^Ws = Pi I ( J 



COS </?i 



This equation can be rewritten as 
(i) 01 d 



(2.31) 



(2.32) 



(2.33) 



Since the string tension in a certain sense characterizes 
the average plaquette value , one may anticipate that this 
quantity is defined (through the duality relations) in the 
terms of the average value of the link £ = -gL- In Zq 

For the 'large' Wilson loop (R m N 2 : L m N 3 ) 



we have 



(2.34) 



N2 



f>2 V3 



6V}A(y) - 



2 cos ip[ 



(2.35) 



and the free energy F = — can be written 



dtp 2 dip 3 
NiJ (2tt) 2 



hi 



(3[RL f dipt cost^i 



N 2 N 3 J 2tt A{(p) 



(2.36) 



The string tension is expressed as 



a^a 2 a 3 



dip 2 dip3 
(2^) 2 



hi 



~ N 2 N 3 J 2tt AMp) 



(2.37) 



It is easy to see that for << 1 we come back to 

1 

get 

<r {1) a 2 a 3 = C(s ) 
where 

r* d V3 



After the integration over Lp\ and 992, one can easily 
C(go-^Q+C(£o+^i) ) (2 38) 



C(*o) 



In (s Q - (3' 3 cos Lp 3 + <J {so - /3' 3 cosip 3 ) 2 - {(3' 2 ) 
In the particular case (f3' n — /?'; a n = a) 



(2.39) 
(2.40) 
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a^a 2 = era 2 



(2.41) 



and we obtain 



, c + O{[3'-(3 ') 2 - s > 3#«f, , 



where the constant Co insignificantly differs for the 'small' 
(co ~ |) and for the 'large' (co ~ |) loops. 
In the more general case 02 = 03 7^ ai 



|h (i "(•>:,>-')■ <^-->> 



It is easy see that in order to obtain the corresponding 
expression for Z (3) the gauge group in given approxima- 



tion one s houl d change only coupling constant in ( 2.45 ) 
given by fl2.12| ) . 

As one may anticipate for 'extremely small' loops 
(R ~ 1; L ~ 1) e.g. sin gle pl aquette, approximate expres- 



sions for 'small' loops (2.31)and (??) become exact 



In a spherical model approximation 'Creutz ratios' are 
the same for 'large' and 'small' Wilson loops 

W (R, L) W (R - 1, L - 1) _ 



W (i?, L - 1) W (R - 1, L) 



^Ws=^/f?V^ (2-44) 



2?r y A{0) ' 

It should be noted, however, that the difference between 
'large' and 'small' Wilson loops has nothing to do with 
finite size effects, it will survive even at an infinite lat- 
tice and reflect alternity between finite Wilson loops and 
infinite ones. At least in given approximation, such dif- 
ference is not so dramatic, as it was forewarned in pl[ . 

Comparison with the MC experiment show that spher- 
ical model predictions qualitively agree with it . Nonana- 
litical, but quite smooth behavior is demonstrated near 
the critical point a (/?'). In the critical region f3' ~ f3' c 
the value cr(P') w a c . In the deep deconfinement region 
/?' << (3' c (where saddle point steadily moved s — > 1 
) a (/?') decrease with g 2 — > smoothly, but too fast if 
compared with the MC experiment. 



D. Spatial string tension 

Here we consider the case of a n = a; ao = a T = a/£ 
and put 

0=-|i Pr = - 2 Z = fc 2 , (2.45) 
9 £ 9 



therefore the expression (2.43) for string tension in spher- 



ical model approximation can be rewritten as 







2 a 



(2.46) 



which in strong coupling region with (2.12)or (2.13) leads 
to 



V<Tfs const ■ , 

a 



(2.47) 



therefore \fa will scale if we demand (3 ~ exp (ax const) 
region the exp 

£e-*; Z(2) 



In weak coupling region the expression ( p. 46} ) with 
(tyjh gives 



(2. 



V^^exp 
and taking into account 



(3 1 — cos 



2tt 

iV 



.9o 



6n In ■ 



1 



a T A L 

obtained in jla] and (l6| we finally get 



with 



2N 

T 



2tt\ , 



l. 



For ^ = 1 the scaling condition e = leads to 
1/2N 



b 



1 — cos 



N , 



N 

—r « .025iV , 

47T^ 



(2.49) 



(2.50) 



(2.51) 



(2.52) 



(2.53) 



which for large enough N agrees with the standard value 
bo = y5 4^t ~ .023^, so all this looks as if the spatial 
string tension on asymmetric lattice (£ << 1) acquires 
'anomalous dimension' e , which disappears when 



III. DISCUSSION 

The application of spherical model in statistical physics 
has long history since the time it has been introduced to 
investigate critical phenomena in the ferromagnet |l7t ] 
and until now ( see, e.g. |23[] , pi]]). Although this model 
is of no direct experimental relevance, it may provide use- 
ful insight since many physical quantities of interest can 
be exactly evaluated with its help. In this context, the 
spherical model is quite a useful tool in providing explicit 
verification of general concepts in critical phenomena, see 

§§§§ 

Recently p3[ it was successfully used for studying 

the transitions between a paramagnetic, a ferromagnetic 
and an ordered incommensurate phase (Lifshitz point). 
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Spherical model approximation helped to find the exact 
scaling function of a system with strongly anisotropic 
scaling. Models of this kind were investigated extensively 
(see recent review in |2j|). 

As it has been established by Stanley (30), there is 
precise correspondence between the spherical model and 
the Heisenberg model. Indeed, consider ad- dimensional 
lattice of N classical spin 

4^/7n*(H^ ( 1 2 ) x 

LI' 



where Uj is a v - dimensional vector of length v s , and 
zjy is corresponding partition function. In the limit 
AT, v — > oo the free energy of this classical Heisenberg 
model is identical to that of the spherical model [j30| . 
Kac and Thompson clarified the situation by proving rig- 
orously that this model gives a surprisingly good result 
for any fixed temperature above and below critical point 
and is independent of the ordering of the limits N — ► oo 
and v — > oo. 

In [p3 it was shown that the spherical model is asso- 
ciated with the matrix model J34J , which has the same 
diagrammatic expansion and saddle points in the pla- 
nar approximation XY model. The analogy with matrix 
models is interesting because it could provide some useful 
technology |34| for solving the XY model. In two recent 
papers ]32|, Parisi et al. introduced and analyzed 
the spherical and XY spin models with frustration to 
test the conjecture that the frustrated deterministic sys- 
tems at low temperature behave similar to some suitably 
chosen spin-glass models with quenched disorder |33| 

Spherical model predicts reasonable values for critical 
exponents Jl^] Moreover, a 'basic' set of exponent rela- 
tions is also satisfied by spherical model for d < 5. 

An advantage of the spherical model is that it satisfac- 
torily describes fluctuations and therefore is suitable for 
the computation of correlation functions either for fixed 
lattice volume or for infinite one. As it was pointed out 
by |35| on lattices with infinite volume the perturbation 
theory (which is in fact a saddle point expansion around 
an ordered state) gives ambiguous results, and certainly 
fails for some border conditions. 



IV. CONCLUSION 

Spatial string tension for 'small' and 'large' loops, as 
one can expect, shows only marginal difference, which 
in Creutz ratio can not be discerned at all. This dif- 
ference, however, is not a lattice artefact and, in princi- 
ple, can be detected in MC experiment. Spatial string 
tension power dependence on lattice spacing, demanded 
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at ( |2.5l| ) , formally does not contradict to condition ob- 
tained by standard renormalisation procedure, ho wever , 
the explicit form of the dependence y/a on £ at (2.51) 
disagrees with it. 
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